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1. Introduction 



Scattering processes involving either initial or final state photons are common in 
electron-positron annihilation, electron-proton collisions and hadron-hadron colli- 
sions. In the initial state, we may be interested in direct processes where the photon 
behaves in a pointlike way, or in processes where the internal structure of the photon 
is probed. Similarly, prompt final state photons may be produced directly at large 
transverse momentum in the hard scattering or from the fragmentation of a large 
transverse momentum jet.^ 

To illustrate the wide range of interesting processes involving photons, let us 
select a few examples relevant to e^e"*", e~p, and hadron-hadron collider experiments 
in the next decade or so. In hadron-hadron collisions, attention is mainly focussed 
on single direct photon production, 

QQ 19, (19 7? 

which is sensitive to the gluonic content of the proton and is an important test of 
perturbative QCD and the pair production of high mass prompt photon pairs, 

qq 77, 99 77 

which is a background to the discovery of a light Higgs boson via its decay into 
photons at the TEVATRON or LHC, 

99^ H ^ 77- 

Important processes in electron-positron collisions include jet photoproduction, 

7? (19, 19 (IQ, 

which is sensitive to the value of the strong coupling constant, as well as the photo- 
production of prompt photons from both the direct 

7? 1(1, 

and resolved processes such as 

9(1 1(1, 

which is sensitive to the gluonic content of the photon. Finally, the proposed gamma- 
gamma option of future linear colliders will, amongst other things, measure dijet 
production, 

77 (IQ 

^Non- prompt photons may also be produced from the decay of a single large transverse momen- 
tum hadron such as the 7r°. 
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which is once again a background to the Higgs boson detection via its decay into 
bottom quarks, 

77 — > if — > bb. 

State-of-the-art theoretical predictions for prompt photon production and the 
vast majority of QED and QCD scattering processes, incorporate corrections at next- 
to-leading order (NLO) in perturbation theory. However, these NLO calculations 
generally exhibit a large sensitivity to the variation of unphysical renormalisation 
and factorization scales as a consequence of truncating the perturbative expansion. 
The inclusion of higher order terms is therefore desirable in order to stabilise the 
theoretical predictions and to reduce the inherent theoretical uncertainties. 

Until recently, a major stumbling block toward a next-to-next-to-leading order 
(NNLO) numerical prediction for scattering processes has been the evaluation of the 
two-loop matrix elements. The first calculation of a two-loop four-point scattering 
amplitude was performed by Bern, Dixon, and Kosower for the maximal-helicity- 
violating gluon-gluon scattering [1]. Subsequently, generic 2^2 scattering matrix 
elements at two-loops have now become tractable for massless particle exchanges 
in the loops and where all of the external particles are on-shell. Smirnov [2] and 
Tausk [3] have provided analytic expansions in e = where D is the dimension, 
in terms of Nielsen polylogarithms for the two most challenging integrals emerg- 
ing; the double-box [2] and the cross-box [3]. At the same time, algorithms based 
on integration by parts [4] and Lorentz invariance [5] recursion relations were also 
developed for the tensor reduction to master integrals of all relevant two-loop topolo- 
gies [6, 7, 8, 9, 10]. 

Already, this technology has been applied to a wide range of physically interesting 
processes. The interference of tree and two-loop graphs (together with the simpler 
self interference of one-loop diagrams) for various processes have now been computed, 
including Bhabha scattering (e+e~ — > e+e^) in the massless electron limit [11] and 
all the QCD 2^2 parton-parton scattering processes {gg — > gg, gg — > qq and 
qq — > qq) [12, 13, 14, 15, 16, 17]. Two-loop hehcity amplitudes have also been derived 
for gluon fusion into photons {gg — > 77) [18], hght-by-hght scattering (77 — > 77) [19] 
and gluon-gluon scattering {gg — > gg) [20] . 

The case where the internal propagators are massless but one external leg is 
off-shell has also been intensively studied, leading to the evaluation of all associated 
planar and non-planar master integrals [21, 22] in terms of two-dimensional harmonic 
polylogarithms [23]. These integrals arise in the decay of an ofi^-shell photon to three 
partons (7* — > qqg) relevant for three jet production in electron-positron annihilation 
and the NNLO matrix elements have been evaluated in Ref . [24] . Even more recently, 
Mellin Barnes integral techniques have been applied to the planar double box diagram 
with four massive and three massless lines, all four legs on the mass shell [25]. 
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In this paper, we present the NNLO virtual corrections for the processes 

q + q ^ 9 + 'y, 
g + g ^ 7 + 7, 
e~ + e+ ^ 7 + 7, 

and the ones related by crossing symmetry or time-reversal. For clarity and cal- 
culational convenience, we decompose the NNLO virtual corrections into the self- 
interference of the one-loop amplitude and the interference of the tree and two-loop 
amplitude. Our results are valid in the limit where the masses of the quarks and 
electrons can be ignored. The Feynman diagrams for the processes in consideration 
are only a subset of those for quark-gluon scattering {qq —>■ gg) , which we presented 
in [15]. However, due to the complicated colour structure of the quark-gluon two- loop 
amplitude and the different flavour content of the processes with photons we choose 
to present them independently. As in [15], we work in dimensional regularisation 
treating all external particles in D dimensions, and renormalise with the MS scheme. 

Our paper is organised as follows. In Section 2 we introduce our notation and 
define the perturbativc expansion of the matrix elements summed over colours and 
spins. In Section 3 we study the singular behavior of the NNLO contributions, and 
verify that it agrees with the general formalism developed by Catani for the infrared 
structure of two-loop amplitudes [26]. For the simpler case of one- loop amplitudes, 
their poles in e can be expressed in terms of the tree amplitude and a colour-charge 
operator I^^^(e), constructed in a universal manner [27]. In the same fashion, the di- 
vergences of the two-loop amplitude can be written as a sum of two terms: the action 
of the l(^)(e) operator on the one-loop amplitude and the action of a new operator 
I*^^^(e) on the tree amplitude. The I*-^''(e) operator includes a renormalisation scheme 
dependent term H^'^'> multiplied by a 1/e pole. Although in [26] it is anticipated 
that if*^^-* can be constructed for any given process in a universal manner, such a 
prescription is not yet available. We give exphcit expressions for I^^^(e) and I^^^(e) 
relevant for each of the processes in Eq. (1.1) vahd in the MS scheme. In Section 4 
we present the finite remainders of the interference of the tree and the two-loop am- 
pfitude after subtraction of the singular poles of Section 3 from the expficit result 
of the two-loop Feynman diagrams. We organize the finite part according to the 
colour and flavour content of the two-loop amplitude. Similarly, in Section 5 we give 
the flnite contributions of the self-interference of the one-loop amplitude. The flnite 
remainders arc expressed in terms of logarithms and polylogarithms which are real 
in the physical domain. Finally, we conclude in Section 6. 
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2. Notation 



We consider the processes 



q{pi) + Qip2) + gips) + lipd 0, 

e~(pi) + e+(p2) + 7(P3) + l{Pi) 0, 
where all particles are incoming with light-like momenta satisfying 



(2.1) 
(2.2) 
(2.3) 



(2.4) 



The amplitudes exhibit infrared and ultraviolet divergences so we work in con- 
ventional dimensional regular isat ion and treat all external particle states in D di- 
mensions. The ultraviolet divergences are renormalised with the MS scheme where 
the bare strong coupling q;° is related to the running coupling as = tts(/x^) at renor- 
malisation scale /i via 



In this expression 



£7 



7 = 0.5772 



Euler constant 



(2.5) 



(2.6) 



is the typical phase-space volume factor in D = 4 — 2e dimensions and /3o, Pi are the 
first two coefficients of the QCD beta function for Np (massless) quark flavours 



lie A - ATrNf 

6 



17C1 - IQCaTrNf - QCfTrNp 
6 



(2.7) 



For an SU {N) gauge theory 



2N 



Ca = N, 



R 



1 

2' 



(2.8) 



In a similar way, the bare QED coupling q;° is related to the running coupling 
a = q;(//'^) at renormalisation scale via 



a 



1 - 



P'o ( ^ 
27r 



+ 



./2 




PS 

2e 



a 
2^ 



+ O {a') 



(2.9) 



where P'q, P[ are now the flrst two coefficients of the QED beta function, 

-2N' 



P'o- 



P'l = -N'p , 



(2.10) 
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and 

f 

with the sum running over the active (massless) fermion flavours. 
The renormahsed amphtudes may be expanded as 



(2.11) 



\M 



OLs_ 

2^ 



\M 



(1) . /tts 

mi) + ^27r 



\Mqq^^) = 4:71 a 



a. 



27r 



a 

2nJ 



(2.12) 



(2.13) 



\M 



Ana 



\M 



(0) 



e~e"'"77 



(1) 

e~e"'"77 



)+(^) lA^l'eO 



+ o 

(214) 

where represents the colour-space vector describing the renormahsed i— loop 

amplitude for the V = qqgj, qqjj, e~e~^jj processes of Eqs. (2.1)- (2.3). The de- 
pendence on the renormalisation scale fi and renormalisation scheme is implicit. 
We denote the squared amplitudes summed over spins and colours by 



{Mqqg^lMqqg.,) = | (? + ? ^ 7 + C/) | ^ = Aqqg^ (s, t, u) , 



(2.15) 
(2.16) 



(Ale-e+77|Ale-e+77) ^ (e e+ ^ 7 7) ^ = A-e+77 («, u) , (2.17) 

where the Mandelstam variables are 



S = (Pl + P2)^ t = {p2+ Psf, 



U 



{Pi+Psf, s + t + u = 0. 



(2.18) 



The squared matrix elements for the crossed processes are obtained by permuting 
the Mandelstam variables and introducing a minus sign for each fermion exchange 
between initial and final states, 

q + q)\'^ ^ Aqqgj{s,t,u), 

Q + 9)\' 

7+?)r 

9 + ^)1^ = -Aqqg^{t,U,s). 



E 1-^(^ + 7 
J2\M{q + ^ 
J2\M{q + g- 
T.\-Mi9 + q- 

E 1-^(7 + 9 



(2.19) 
(2.20) 
(2.21) 
(2.22) 

(2.23) 



Similarly, 



E 1-^(7 + 7 
El-M(? + 7- 
E 1-^(7 + 



■g+g)r 

? + 7)r 

7+?)r 



(2.24) 
(2.25) 
(2.26) 
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and 



A. 



[s,t,u), 



-A. 



'u,t,s). 



(2.27) 

(2.28) 
(2.29) 



The functions A-p{s,t,u) are symmetric under the exchange of t and u and can 
be expanded perturbatively to yield, 



Aqqg-y {s, t, u) = IGii^aas 



+ (|)'-4Sr(M,«) + 0(aX) 



(2.30) 



-4j°,^i,«) + (^j<^;?(s,*,t<) 



+ (|)'A"^f(..t.«) + o(„,V) 



and 



(s, t, m) = IQiT^a^ 



■e~e"'"77 



■e-e+77 ^ ' ' l27r/ ^~^''''>")' 



+ (^)'A™.';°,(M,«) + o(a') 



At leading-order (LO) the self-interference of the tree amplitudes is 

A'f{s,t,u) = {Mf\Mf), 
so that for each process we have 

A];g.^{s,t,u) = NCFr{s,t,u), 

(s,t,M) = Nr{s,t,u), 



9977 

) 

••e-e+77 

where we have defined the tree-type structure 



^/ V ^ ^ I u t s 

T(s,t,u) = 8(1 - e) - H e— 

^ ^ ^ ^ \ t u tu 



(2.31) 



(2.32) 



(2.33) 



(2.34) 
(2.35) 
(2.36) 



(2.37) 



The next-to-leading order (NLO) term consists of the interference of the tree 
and the one-loop amplitudes 



(0)1 A^(l)\ 



(2.38) 
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and their expansion in e up to (9 (e) was given in Ref. [28]. As will be discussed in 
Section 3, the singular structure of the two-loop amplitude is expressed in terms of 
the tree amplitude multiplied by a singular operator of 0{l/e*) and the one-loop 
amplitude multiplied by another operator diverging as Therefore, one needs to 
know the expansion of the one-loop amplitude up to and including the O (e^) term. 
In Section 3 we give expressions for {M^^\Ai^^) vahd in all kinematic regions and to 
all orders in e, in terms of two one-loop master integrals. The e-expansions of these 
master integrals to the appropriate order are given in the Appendix C. 

At next-to-next-to-leading-order (NNLO) contributions from the self-interference 
of the one-loop amplitude and the interference of the tree and the two-loop amplitude 
must be taken into account, so that 

A^'^'^^'is, t, U) = ^NNLO(lxl) ^ ^J^LO(OX2)^^^ ^2.39) 



with 



and 



AT''^''''\s,t,u)^{M^i^\M^i^), (2.40) 



^NNLO(ox2)^^^^^^^ = {M^^W?) + {M^?\M^P). (2.41) 

In the following sections, we present expressions for the infrared singular and fi- 
nite contributions to A^^^'~'. In Section 3 we write down the singular parts ac- 
cording to the formalism of Catani [26]. The pure two- loop finite contributions 

are described in Sec. 4 and the finite remainders of the self- 
interference of the one-loop amplitude A^^^'^^^^^\s,t,u) are described in Sec. 5. 

As in Refs. [12, 13, 14, 15, 16, 17], we use QGRAF [29] to produce the tree, one 
and two- loop Feynman diagrams to construct \M.^). We then project by or 
( Alp^ [ and perform the summation over colours and spins. It should be noted that 
when summing over the gluon polarisation, we ensure that the polarisation states 
are transversal (i.e. physical) by using an axial gauge 

Ee^C = -^- + ^^^^^^^^, (2.42) 

where n is an arbitrary light-like 4- vector. For simplicity, we choose — P4- Finally, 
the trace over the Dirac matrices is carried out in D dimensions using conventional 
dimensional regularisation. It is then straightforward to identify the scalar and tensor 
integrals present and replace them with combinations of the basis set of master 
integrals using the tensor reduction of two-loop integrals described in [6, 7, 8], based 
on integration-by-parts [4] and Lorentz invariance [5] identities. The final result is a 
combination of master integrals in D — A — 2e for which the expansions around e = 
are given in [2, 3, 6, 7, 8, 9, 10, 30, 31]. 
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3. Infrared Pole Structure 



We further decompose the one-loop self-interference and the two-loop contributions 
as a sum of singular and finite terms, 



^NNLO(lxl) . 



s, t, u) - pr., t, u) + t, u) 



(3.1) 



and 



(3.2) 

for each of the processes V — qqg'j, 9977, e~e"'"77. V^j^g^^ and VJI^J^^ contain infrared 
singularities that will be analytically canceled by the infrared singularities occurring 
in radiative processes of the same order (ultraviolet divergences having already being 
removed by renormalisation) . J-lnUe,-p '^nUe,'p remainders which are finite 

as e ^ 0. 

The poles of the one-loop amplitude self-interference can be written in terms 
of a universal operator I*-^''(e) acting on the colour-space of the amplitude. Due to 
the simple colour structure of the processes we study, the action of I*^^'*(e) factorises 
yielding. 



)\ {MP\MP) + 2Re{l^i\e)HM^^Wv^)} (3-3) 



where 



id) 



on 



2T{1 - e) 




3 






2Ne 












(3.4) 



and 




(3.6) 

The generic form of the I*^^^ (e) operator was found by Catani and Seymour [27] and it 
was derived for the general one-loop QCD amplitude by integrating the real radiation 
graphs of the same order in perturbation series in the one-particle unresolved limit. 
For the processes qqg'y and 9977 we consider only NNLO QCD corrections, therefore 
in Eq. (3.5) there are no contributions from soft or coUinear photon emission at the 
same order in a^. Note that on expanding the I^^^(e) operator, imaginary parts are 
generated, the sign of which are fixed by the small imaginary part +iO assigned 
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to each Mandelstam variable. Terms that involve the hermitian conjugate of this 
operator are to be modified accordingly. 

The renormalised interference of the tree and the one-loop amplitudes can be 
expressed in terms of the one-loop bubble graph (Bub) and the one-loop box integral 
in D = 6 — 2e dimensions (Box^), as follows 



C 

Cpfi {s, t,u)-^ {/i (s, t, u) + /2 (s, t, u)} 



{MZ,\M\',>,,) = NCf 

{M%^\M%^) = NCfMs, t, u) + {t^ u) (3.8) 



- Ye^Mfl,\M^Sj, (3.7) 



and 



(•^i-U7l-^i-U7) = ^' + ^) - 7 (-^l-UTl-^e-UT)' (3-9) 

where we have used 

, ^ _ 8 (1 - 2e) [g^ + -eis"^ + + u'^) + {t^ + - st) - e^st] ^ ^ 
Ji[s,t,u) — ^ rSox [s,t) 

eut 

4(l-.)[3. + t + .(5t + 2«)^BubW (3.10) 



t 

and 

f2{s,t,u) = ^ ^-^^ ^ ^ ^Box6(i,«) 

4 (1 - e)' [2 (^2 + «2) _ ^ _ _ ^2^^ 



Bub(t). (3.11) 

These expressions are valid in all kinematic regions and for all orders in e. However, 
to evaluate them in a particular region, the one- loop master integrals Bub and Box^, 
must be expanded as a series in e (see Appendix C). 

The infrared poles of the interference of the tree and the two-loop amphtudes 
follow a generic formula developed by Catani [26], such that 

VZ%{s,t,u) = 2Re{l«(e)(Mf |M«) +lg)(e)(Mf |a4°))} (3.12) 

where once again the simple colour structure of the processes under consideration 
allows the action of I*-^''(e) and I*-^'*(e) to be factrorised. For the QCD processes 
^ = QQ91y we have 

lS?)(e) = -il<;)(e) (lW(e) + + tl^^ + k) I«(2e) 

H!^> (3.13) 



4er(l - e) 
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with 

and for the QED process V — e~e"'"77 we have 

^i-U. (3-15) 



4er(l - e) ^ ^^^T 
with 

X' = -^N'^. (3.16) 

The process and renormahsation scheme dependent H^'^^ constants are related 
to the colour space operator H*^^^ The full structure of H*^^^ is not known at present, 
although it certainly contains non-trivial colour structures [20] that have been in- 
vestigated in the case of gluon-gluon scattering. However, because of the projection 
by the tree-level amplitude and the summation over colours we are only sensitive to 
the trivial colour part of H*-^-*, and wc find that the constants H^'^'' can be written 
as a sum of factors associated with each of the external legs. We derive them di- 
rectly from the 1/e term of the explicit expansion in e of the two- loop amplitudes. 
In conventional dimensional regularisation and the MS scheme, we find that for the 
processes under consideration 

i/g), = iff) + i7f (3.17) 

^S-77 = ^f + ^ (3-18) 
Hil,^^^HS} + HiU2H(^\ (3.19) 

where for an external quark or an antiquark we find 

q q ^2 8y ^ V 2 ^ 216 48 / 

+ (-| + ^) TnN^Cp (3.20) 
and for each external gluon 

These are the same factors as found in [26, 12, 13, 15, 16, 24]. Similarly, for each 
external electron or positron we find 
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while for the external photon 

H^''> = |iV;^ + iV;. (3.23) 

It should be noted that H^^^ and H^^^ can be obtained from the expressions for 
Hf\Hf^ and H'p accordingly, by taking the limit Tr ^ 1, Cf ^ 1 and Ca ^ 0. 
It should also be emphasized that contributions of O (e) to H^'^^ are undetermined 
at present. 



4. Finite two-loop contributions 

In this section, wc give explicit expressions for the finite remainder of the two-loop 
contributions J^inul^-p defined as, 

^n^%(«, u) = t, u) - VSris, t, u). (4.1) 

Note that the V^^i^/^i^s, t, u) is expanded through to O (1) and therefore contains finite 
as well as singular contributions. Using the standard polylogarithm identities [38], 
we express our results in terms of a basis set of logarithms and polylogarithms^ with 
arguments x, 1 — x and {x — l)/x, where 

t u u X — 1 ^. 
x = — , y= — = l-x, z^--^ . (4.3) 

S S t X 

In the physical region s > and t, m < 0, our basis set of functions are all real. For 
convenience, we also introduce the following logarithms 

X = log y = log 

(4.4) 

where is the rcnormalisation scale. 

Prom Eqs. (2.15)-(2.29) we see that we need to evaluate the function T^ftl,-p '^i^^^ 
arguments (s,i,w), {u,t,s) and {t,u,s). Since t and u are both negative they can 
be interchanged leaving the polylogarithms and logarithms of our basis well defined 
and real. Therefore the crossing symmetry {t ^ u) can be performed trivially. It is a 
more involved operation obtain expressions valid under the exchange of s with either 
t or since the logarithms and polylogarithms may acquire imaginary parts and 

^As usual, the polylogarithms Li„(w) are defined by 

/■"' dt 

Li„(u;) = / — Li„_i(t) forn = 2,3,4 
Jo ^ 

U2{w) = - ^log(l-i). (4.2) 
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(a) 



> — t/VW^^ — tAW 



AAAAA 




> — t/WV 



WW (b) 



< — vvw 



Figure 1: Feynman diagrams representative of (a) "light-by- light" scattering and (b) 
"abelian gluon" exchange. These graphs contribute to Ac and Be respectively. Depending 
on the process, the zigzag lines represent either photons or gluons. 



appropriate analytic continuations need to be considered. We therefore give explicit 
expressions for the "s-channel" function J^°j^g ,p(s, t, m) and the "w-channel" function 



^0x2 
•'inite,'p 



(u,t,s) that are directly valid in the physical region s > and u.t < 0. 
The "t-channel", J-^^^^'^ ^{t,u, s) can be obtained from the -u— channel expression by 
exchanging t and u. Note that here we define a "channel" according to the first 
argument of the function J-inite,^,- 

We organize our results according to the colour and flavour structure of the 
two-loop amplitude, so that for the three different processes, we have 



J=l^^'^ _ = 2NCf 

t'li't'i'f^^qqg'y- q 



+ClB, + C\C, + CfCaD^,, + NfCfE^,, + NfCaE^.,^ + NIF^,, 



^0x2 
"^inite, 



2N 



(JIJ77; c 



^ Ql TrCfA, + ClB, + CfCaD2;c + NfCfE^,, 



(4.5) 
(4.6) 



and 



^0x2 
mite, c+ ^ 



(4.7) 



where the subscript c = s, -u denotes the channel. refers to the charge of each 
of the q active masslcss quark flavours for the QCD processes and cquivalently Qf 
denotes the charge of each of the / active massless fcrmion flavours participating 
in the e~e"'"77 process. We see that the "light-by-light" contribution, A^, due to 
internal fermion box graphs and the "abelian gluon" contribution, i?c, illustrated in 
Fig. 1, occur in all three processes. 

The values of Ac, Be, Cc, Di.c, E^c and F^c are presented in Appendix A.l for 
the s-channel and in Appendix A. 2 for the li-channel. 
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5. Finite one-loop contributions 

In this section, we give explicit expressions for the finite remainder of the one-loop 
self-interference contributions ^^j^e,-p defined as, 

(5.1) 

Note again that the 'Pji^g^^{s,t,u) is expanded through to O {!) and therefore con- 
tains finite as well as singular contributions. In fact, as in [15, 14, 17], all polyloga- 
rithms are collected in the singular part ('PoLs,^p)> leaving finite remainders with only 
logarithms and constants such as tt^ and Cs- This is because in the expansion of the 
box integral in D = 6 dimensions, (see Appendix C) the polylogarithms appear in 
the O (e) and O (e^) terms. In order for these polylogarithms to contribute at O (1), 
they must be multiplied by an infrared singular term and are therefore contained 
in Pjes^-p- The finite part of the one-loop self-interference contribution due to the 
interference of one box graph with another only collects the purely logarithmic O (1) 
terms in each. 

Specifically we find, 

CpGi-c + CaCfG2;c + C^G^-c 

+NfCfX,., + NfCaX2;c + A^|^3;c 

^initiqq^^-c^ ■^^'^Gl-c (5-3) 

and 

^nule-e^,,^ = ^^'^c + N'^X,,, + N'^' X,,,, (5.4) 

where the subscript c — s,u denotes the channel. 

The values of G^c and Xj-c are presented in Appendix B.l for the s-channel and 
in Appendix B.2 for the zi-channel. 



(5.2) 
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6. Conclusions 



In this paper we calculated the NNLO QED and QCD virtual corrections for a 
range of 2 — > 2 massless scattering processes with external state containing photons, 
namely qqg^i qq'y^, e~e^^^ and the processes related by crossing symmetry and time 
reversal, in the high energy limit where the fermion masses can safely be ignored. 
We used conventional dimensional regularisation and to remove the UV divergences 
we renormalised using the MS scheme. The renormalised amplitudes are infrared 
divergent and contain poles to O (1/e^). In fact we were able to check our results by 
comparing the infrared structure of the calculated amplitudes with the prediction of 
Catani's formalism [26] for the infrared structure of generic one- and two-loop QCD 
amplitudes. Catani's method does not determine the 1/e poles exactly, but suggests 
that the undetermined non-logarithmic contribution if^^^ can be extracted from a few 
basic two-loop amplitudes. Prom the amplitudes calculated in this paper together 
with Refs. [15, 16, 24] we find that, when summing over colours, each external-leg 
contributes independently to the leading order term in e of H'^'^\ The contributions 
of quarks, anti-quarks and gluons and their QED analogues are given in Section 3 
for the MS scheme. 

The main result of our paper is the finite remainders of the NNLO virtual correc- 
tions after we subtract the pole structures of Catani's formalism expanded through 
to C (1). In Section 4 we gave the finite remainders of the interference of the tree 
and the two-loop amplitude and in Section 5 we presented the finite remainders of 
the self-interference of the one-loop amphtude for each of the processes under con- 
sideration. The results are expressed in terms of poly logarithms and logarithms that 
are real in the physical region. The renormalised amplitudes for the three processes 
WQl: QQll cind e^e+77, difi^cr with respect to group invariants, the renormalisation 
procedure, and their flavour content. 

The aim, of course, are more precise predictions for the basic scattering pro- 
cesses. Initial studies suggest that at NNLO an accuracy of a few percent is attain- 
able for strong interaction processes. However, much work is needed to accomplish 
this goal. An important ingredient is a systematic procedure for analytically car- 
rying through the cancellation of the infrared singularities present in the virtual 
contributions against the contributions from the one-loop 2 — > 3 processes when one 
particle is unresolved and the trcc-lcvcl 2-^4 processes when two particles arc unre- 
solved. Such a method has not yet been established, although the single and double 
unresolved limits of the relevant matrix elements are well known [32, 33, 34, 35]. 
In fact, many of the analytic phase space integrations for the double unresolved 
and single unresolved loop contributions have already been studied in the context of 
e+e~ — > photon -|- jet at C(qQs) [36] and Higgs production in hadron colliders [37]. A 
further complication may be due to initial state radiation where the three loop split- 
ting functions [39, 40, 41, 42, 43] are needed to extract parton density functions [44] 
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at an accuracy matching that of the hard scattering matrix element. Nevertheless, 
it is an important task to achieve more reliable theoretical calculations that can 
take advantage of the improving experimental data and make a better test of the 
underlying physics at short distances. 

Acknowledgements 

C.A. thanks the theory group at Fermilab for their kind hospitahty during the early 
stages of this work. C.A. would like to thank Zvi Bern and Lance Dixon for useful 
conversations. We gratefully acknowledge the support of the British Council and 
German Academic Exchange Service under ARC project 1050. This work was sup- 
ported in part by the EU Fourth Framework Programme 'Training and Mobility of 
Researchers', Network 'Quantum Chromodynamics and the Deep Structure of Ele- 
mentary Particles', contract FMRX-CT98-0194 (DG-12-MIHT) and by the National 
Science Foundation, grant PHY9722101. 



15 



A. Finite two-loop contributions 

A.l s-channel 

In this section we list the finite contributions for one-and two-loop contributions in 
the s-channel as defined in Eqs. (4.5), (4.6), (4.7), (5.2), (5.3) and (5.4). 
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In this section we list the finite contributions for one-and two-loop contributions in 
the li-channel as defined in Eqs. (4.5), (4.6), (4.7), (5.2), (5.3) and (5.4). 
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B. Finite one-loop contributions 
B.l s-channel 



Gi-s 
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B.2 ii-channel 
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C. One- loop master integrals 



In this appendix, we list the expansions for the one-loop box integrals in D — 6 — 2e. 
We remain in the physical region s > 0, u,t < 0, and write coefficients in terms of 
logarithms and polylogarithms that are real in this domain. More precisely, we use 
the notation of Eqs. (4.3) and (4.4) to define the arguments of the logarithms and 
polylogarithms. The polylogarithms are defined as in Eq. (4.2). 
We find that the box integrals have the expansion. 



Box^ {u,t) 



e^^r(l + e)r(l-6)^ 
2sr(l-2e)(l-2e) V s 

+2e 



{X - yy + 

2 



-2e' 



Li3(a:)-XLi2(x)--X ^ 

Li4(a;) + rLi3(a;) - lx^U2{x) - - ^X^Y + ]x^Y^ 
Z o D 4 



TT^ TT^ TT^' 

.!LX^ - -XY - — 

4 3 45 



+ {u^t)\ + 0{e'), (C.l) 



and 



(-2Li3(x) + 2XU2{x) - + 2YX^ - n'^X + 2(3) 



TT 



+i7r 2Li2(a;) + iYX - X' - — 



2U4iz) + 2Lu{y) - 2FLi3(a;) - 2XLk{y) + {2XY - X^ - 7r^)U2{x) 
+\x^ - ^X^Y + ^X^Y^ + '^TT^X^ - 2n^XY + 2FC3 + ^ttM 

+i7r(^-2U3{x) - 2Li3(y) + 2YU2{x) + ^X^ - 2X^Y + 3XY^ 



-Y + 2Cs 



+ O (f) . 



(C.2) 



Box^(s,-u) is obtained from Eq. (C.2) by exchanging u and t. 

Finally, the one-loop bubble integral in D = 4 — 2e dimensions is given by 



Bub(s) 



e--r(i + r)r(i-r)V 

r (2 - 2e) e s 



(C.3) 
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